A study of the correlations between jet quenching observables at RHIC 
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Focusing on four types of correlation plots, Raa vs. V2, Raa vs. 7aa, Iaa vs. ^2** and V2 vs. 
iij** , we demonstrate how the centrality dependence of correlations between multiple jet quenching 
observables provide valuable insight into the energy loss mechanism in a quark-gluon plasma. In 
particular we find that a qualitative energy loss model gives a good description of Raa vs. V2 only 
when we take AE l^ and a medium geometry generated by a model of the Color Glass Condensate. 
This same AE ~ l^ model also qualitatively describes the trigger pr dependence of Raa vs. Iaa 
data and makes novel predictions for the centrality dependence for this 7?aa vs. Jaa correlation. 
Current data suggests, albeit with extremely large uncertainty, that ^2** ^ V2, a correlation that 
is difHcult to reproduce in current energy loss models. 
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Introduction: The combination of theoretical pre- 
dictions and experimental measurements of high-p^ jet 
observables provides a unique basis for determining the 
properties of the strongly-interacting quark gluon plasma 
(sQGP) created in Au-|-Au collisions at the Relativistic 
Heavy Ion Collider (RHIC) [1-3]. After nearly a decade 
long effort, jet quenching via final state partonic interac- 
tions, as an experimental phenomenon, has been firmly 
established at RHIC [3]. The challenge for theory is to 
find an energy loss model built on first principles deriva- 
tions that simultaneously describes the known observ- 
ables. Currently no such model exists, and there is a de- 
bate over the exact, or even dominant, energy loss mech- 
anism at work in the sQGP (sec, e.g., [3-5]). 

Four observables of interest in leading particle quench- 
ing physics are single hadron suppression Raa, V2 (one 
half the coefficient of the cos(2(/)s) term in the Fourier ex- 
pansion of the suppression relative to the reaction plane 
(RP) RAA{<f>s = </> — *Rp)), di- hadron suppression Iaa, 
and . These observables are interesting because they 
probe the same energy loss processes in a heavy ion colli- 
sion, but with different underlying parton spectra and/or 
path length 'T' dependencies. For example Raa at dif- 
ferent (ps has identical input parton spectra but explores 
different path lengths; on the other hand, Iaa probes a 
harder input spectrum and a different set of paths com- 
pared to Raa- The importance of using multiple observ- 
ables to constrain the possible energy loss mechanism 
in heavy ion collisions is well known [G-10]. However, 
other than one previous publication which examined the 
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centrality dependence of Raa vs V2 [11], the comparison 
between theory and data was made one observable at a 
time and usually as a function of px for one centrality 
selection. 

In this work we argue that correlating these observ- 
ables directly against each other and studying the cen- 
trality dependence of the correlation provides novel in- 
sights into the high-p^^ energy loss mechanism. We pro- 
pose four types of correlations that can be studied exper- 
imentally at RHIC and the LHC: Raa vs. V2, Raa vs. 
Iaa, Iaa vs. v^^'^ and V2 vs. W2**- We shall first give a 
brief overview of the experimental measurements of each 
observable at RHIC. We then explore the main features 
of the four correlations revealed from the experimental 
data. Using a jet absorption model, we demonstrate the 
importance of these correlations by exploring their sen- 
sitivities on the geometry and parton spectra shape. We 
conclude with a discussion of how these correlations can 
be used to disentangle the dependence of energy loss 
from the collision geometry and parton spectra. 

Overview of experimental results and theoret- 
ical comparisons: The most precise measurements of 
high-pT single hadron suppression and anisotropy were 
carried out by the PHENIX experiment using 7r° mesons 
[12, 13], reaching pT ~ 20 GeV/c for Raa and beyond 
10 GeV/c for V2- Raa is defined as 

r'' ('^ ^ M — d^PTiPT;<Ps,b) 
I^AAiPT, (ps, b) = — rfjVPP^i>+x , (1) 

where Ni,in{b) is the number of binary (hard scatter- 
ing pp-\ikc) collisions at impact parameter b, and V2 = 
/ d(f>s Raa{<I>s)cos{2(I)s)/ j d(l)sRAA{4>s)- The i?AA shows 
an almost pT-independent factor of 5 suppression in cen- 
tral collisions for pT > 4 GeV/c. The V2 drops from 3 
to 7 GeV/c, but remains positive at higher pT- Current 
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jet quenching models based on the pQCD framework, 
when tuned to Raa data, significantly under-predict the 
172 [13]. In contrast, non-perturbative approaches, for 
example those based on AdS/CFT gauge gravity dual- 
ity [14], seem to work well. The data seem to prefer 
the AE path length dependence, a result based 

on AdS/CFT [7], as opposed to the quadratic depen- 
dence /S.E ^ P predicted radiative energy loss predicted 
by pQCD [15]. Alternatively, a simultaneous descrip- 
tion of Raa and V2 may also be achieved via a late-stage 
non-perturbative effect near the QCD confinement tran- 
sition [11, 16]. 

The suppression of the away-side jet is quantified by 
Iaa , the ratio of the per-trigger yield (away-side jet mul- 
tiplicity normalized by number of triggers) in Au-|-Au 
collisions to that in p-l-p collisions. Pure geometrical 
considerations would imply Iaa < Raa, due to a longer 
path length traversed by the away-side jet. But recent 
PHENIX [17] and STAR [18] measurements show that 
Iaa is constant for associated hadron > 3 GeV/c, 
within the current experimental uncertainties, and this 
constant level is above the Raa for the trigger hadrons, 
i.e. Iaa > Raa (see Fig. 1). Furthermore, the constant 
level of Iaa increases for higher trigger p^. This result 
can be qualitatively explained by the bias of the away- 
side jet energy by the trigger pT'- as we show with a 
PYTHIA simulation in Fig. 2, the initial away-side jet 
spectra becomes harder as higher pT triggers are re- 
quired. As we discuss in more detail below, the harder 
the input spectrum, the larger the fractional energy loss 
is required for the same Iaa value. The ACHNS model 
[19] results, constrained by the Raa data, are incom- 
patible with the Iaa values shown in Fig. 1; while the 
ZOWW model [20] seems to describe the Iaa data alone 
shown in Fig. 1, it too fails at simultaneously describing 
both Raa and Iaa [21]. 

PHENIX [22] recently reported the first measure- 
ment of the anisotropy of away-side suppression, -y2** = 
J d4>slAA{4>s)cos{2(j)s)/ J d(j)slAA{(l>s)- The away-side 
yield shows a strong variation with angle of the trigger 
relative to the RP. This variation is much larger than 
that for inclusive tt*^ in the same trigger p^ range. The 
current measurement is statistics limited; however the re- 
sult is tantalizing as energy loss models usually predict 
much smaller Wj** [22]. 

The jet absorption model: We use a model 
from [23, 24] to investigate correlations between the four 
observables, and to check the sensitivities of these corre- 
lations to the collision geometry and I dependence of the 
energy loss. The model is based on a nai've jet absorp- 
tion picture where the fractional energy loss of a high px 
particle is proportional to a line integral / through the 
medium, e = kl, where e = 1 — p^/pl^. For a power law 
production spectrum with index of n, dN/(PpT P^"' 
i?AA is related to the fractional energy loss e = AE/E 
via [25, 26]: 

Raa = ((1 - e)"-') « {e'^'"^^') « (e"'^^), (2) 
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FIG. 1: (Color online) (Top panels) PHENIX tv° - h corre- 
lation results with 7 < Pr < 12 GeV/c and 0.5 < Pr < 7 
GeV/c [17]. (Bottom panel) the STAR jet-h correlation 
results with reconstructed trigger jet momentum in 10 < 
PtXc < 40 GeV/c [18]. Both are for the 0-20% Au-^Au 
centrality bin. 
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FIG. 2: (Color online) PYTHIA simulation of the away-side 
per-trigger charged hadron yield spectra in PHENIX rj ac- 
ceptance for three different charged hadron trigger pr ranges. 
The away-side yields are parameterized with a power law func- 
tion. 

where k = (n — 2)k, and (• • • ) indicates an average 
over the binary collision profile. The line integral / is 
calculated as /i = J p dl or I2 = f p Idl. The for- 
mer corresponds to a quadratic dependence of energy 
loss {dE ^ Idl) in a longitudinally expanding medium 
(p(r) ^ 1/t ^ l/l), while the latter corresponds to a cu- 
bic dependence {dE ^ Pdl) of energy loss in a longitudi- 
nally expanding medium. Up to slowly varying logarith- 
mic factors the interference of the unmodified vacuum 
radiation associated with the production of a high-py 
parton with the medium-induced bremsstrahlung radia- 
tion in the deep Landau-Pomeranchuk-Migdal region — 
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which one expects with the ordering of length scales 
<^ Xjnfp ^ Hn a weakly-coupled QGP described by 
Hard Thermal Loop pQCD weakling interacting with the 
high-pT parent parton that we might expect in Au+Au 
collisions at RHIC — yields a fractional energy loss that 
scales with the square of the pathlength, de oc Idl. For the 
not-too-large fractional energy losses at RHIC, (e) ~ 0.2 
for Raa ~ 0.2, the exponential absorption model is a 
reasonable approximation to the 1 — e Jacobian expect 
for pQCD energy loss. On the other hand, under the 
assumption that all the couplings between the sQGP 
and the high-p^ parton are very large and the dominant 
physics can be well approximated using the AdS/CFT 
conjecture, then the thermalization distance for a light 
high-pT probe parton scales as E'^/'^. The exponential 
model can again be used, in this case capturing the 
physics of the probability of escape for the strongly cou- 
pled high-pT particles. We model the medium density p 
either by the participant density profile from Glauber ge- 
ometry or gluon density profile from CGC geometry [27]. 
The dominant effect of event- by- event fluctuations in the 
sQGP arc included in a medium rotation procedure [28]. 
Comparison of the jet absorption model to data is a rea- 
sonable first step when examining the physics of the cen- 
trality dependence of the correlations investigated in this 
paper as it captures both the pathlength dependence and 
medium geometry effects. 

K is the only free parameter in our energy loss model; 
we tune it to reproduce Raa ^ 0.18 for 0-5% most cen- 
tral TT*^. Once K is fixed, we then predict the central- 
ity dependence of Raa in 5% centrality increments, as 
weU as Iaa, and Wj**- The k values for the four 
cases (the combinations of l^, P and Glauber, CGC me- 
dia) are summarized in Table I. Note that for a given 
"I" dependence, the suppression level is essentially con- 
trolled by the product of k and the average matter den- 
sity (pmodium) = / p(x)^d^a;/ / p{K)d'^x in the 0-5% cen- 
trality bin. In general the «;(pmcdium) for CGC geometry 
is slightly larger than Glauber geometry, primarily be- 
cause the former has a smaller matter profile [24], while 
both geometries are assumed to have the same binary 
collision profile. 



TABLE I: K, average matter density (pmcdium) = / p{5i)^d^x/ 
J p{x)d^x, and the product of the two in the 0%-5% Au+Au 
centrality bin, for the four cases calculated in our study. 
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If we assume that the di-hadron production spectrum 
is also a power law, dN / {d^p'^d^p'i^) ~ (Pt)~"° (Pt)~"' 



then 

Iaa = ((1 - ea)""-'(l - et)"'-')/i?AA 

« (e-'^™-^"e-«^*)/^AA. (3) 

Since Kocrit — 2 = n — 2 according to Eq. 2, the effec- 
tive Kaway for away-sidc jets should be smaller due to a 
smaller n^. Fig. 2 shows that increasing the momentum 
of the trigger particle from 4 — 5 GeV/c to 9 — 12 GeV/c 
yields a reduction in the power law for the away-side 
spectrum from — 5.1 to = 3.6. One may effec- 
tively model this reduction in away-side input spectrum 
in our absorption model by approximating 

^away ^ ^5 (4) 

hence the effective strength of the energy loss is reduced 
by a factor between 2 and 4 for the trigger particle mo- 
mentum ranges currently measured. One may readily 
see that as the trigger momentum range is increased, the 
effects of energy loss on the suppression of particles is 
reduced; Iaa is not simply smaller than Raa due to the 
longer pathlengths that result from the trigger bias in the 
di-hadron measurement. 

Results: Fig. 3 (a) shows the predicted correlation 
between i?AA vs. V2 from the jet absorption model over 
the full centrality range. The calculations appear to show 
little dependence on the assumed geometry (more later), 
but clearly V2 increases dramatically from P to depen- 
dence. The dependence agrees with data well, implying 
that it can simultaneously describe both i?AA and U2, a 
conclusion already made in [24]. 

We know that the low-p-r V2 is observed to scale with 
eccentricity (e) [29]. It was shown in [24] that the jet 
quenching V2 also approximately scales with e. Thus it 
is instructive to plot i?AA versus the reduced quantity 
V2.r ~ as shown in Fig. 3 (b). The data now appear 
as two sets of points, filled black circles and open red cir- 
cles, corresponding to Glauber geometry or CGC geome- 
try, respectively. They both indicate an anti-correlation 
with i?AA, that is a large V2,r corresponds to a small 
i?AA and vice versa. Similar trends are also shown by 
the calculations: as quenching becomes stronger the sur- 
viving jets further amplify the initial asymmetry. In 
Fig. 3 (c) we show the centrality binned theoretical re- 
sults as open black diamonds and open red crosses for the 

Glauber and CGC medium results respectively; these 
theory points should be directly compared to the filled 
black circle data points. 

Note that while the e ^ models with either a CGC 
medium or Glauber medium appear to describe the i?AA 
vs. V2 data in Fig. 3 (a) well, only the cubic model with 
the CGC medium describes the Raa vs. V2,r data shown 
in Fig. 3 (b). As shown in Fig. 3 (c), this is because 
the f' model with Glauber medium does not describe the 
data at the correct centrality bin whereas the model 
with CGC medium does. Eccentricity is a centrality de- 
pendent quantity (as discussed in detail in [24] the CGC 



4 



" /" Glauber I Glauber 

I' CGC l' CGC 

71° Vj 0-50% 
in 10% steps 
p e[6-91 GeV/c 



< 
DC 



0.5- 




(a) 



. '' "i/'^cc.c 



-I' v,/e 




(b) 



tC 0.4 



0.2 



0.05 0.1 0.15 ) 



I' Glauber 

* CGC 




0.05 



0.1 



FIG. 3: (Color online) (a) Centrality dependence of -Raa vs. 
V2 from the jet absorption model in 5% centrality steps and 
data (solid circles); only statistical errors are shown, (b) The 
same data and calculations, except V2 is divided by the ec- 
centricity, (c) The same as in (a) but zoomed in and with 
centrality binning explicitly shown for the theoretical results. 



geometry is smaller relative to the Glauber geometry), 
and normalizing with respect to eccentricity has empha- 
sized this mismatch in the centrality-binned results for 
theory and data. Thus i?AA versus U2/e can better illus- 
trate the relation between jet quenching and azimuthal 
anisotropy, indicating here that our model only describes 
the -Raa vs. V2 data for AdS/CFT-like energy loss in a 
CGC medium. 

Fig. 4 shows the correlation between i?AA and /aa 
from the jet absorption model, calculated over the full 
centrality range. The model results are compared with 
STAR and PHENIX data from Fig. 1 (0-20% central- 
ity bin) with their /aa values integrated over > 3 
GeV/c (where /aa is flat). In Fig. 4 (a), the data points 
have roughly the same /?aa value, but are spread out 
in /aa for different trigger momenta, p^. The reason, 
as explained in the discussion of Fig. 1 & Fig. 2, is that 
/aa depends not only on the path length but also on 
the shape of the input spectra, and the larger the trig- 
ger momentum ptp the harder the away-side spectrum. 
Fig. 4 (a) also shows that /aa < -Raa when only the 
path length effect is included (i.e. we take Kaway = k). 
We then attempt to model the effect of the trigger bias 
on the hardening of the away-side spectrum (see Fig. 2) 
by using Eq. (4), which yields Kaway = k/2 in Fig. 4 



(b) (p^ G [4-5] GeV/c), 
(p^ e [7-9] GeV/c), and 



k/3 in Fig. 4 (c) 
k/4 in Fig. 4 (d) 
(py G [9 — 12]). The toy model improves its agreement 
with the PHENIX data when both the path length and 
spectral dependencies are included. 

In particular the AdS/CFT-like energy loss model 
that described the /?aa vs. V2 data so well appears to de- 
scribe the /?AA vs. /aa data to within about 2-3 standard 
deviations. One also again sees that the CGC medium 
yields results whose centrality dependence is in slightly 
better agreement with the data than the results from 
the Glauber medium. However it is clear that the P 
models systematically under-predict the /aa data. On 
the other hand the P models tend to disagree more on 
the level of 1 standard deviation. That the models 
tend to describe the normalization and correlation-but 
not anisotropy-of /?aa and /aa might suggest the im- 
portance of hadronization or flow-coupling effects that 
are neglected in our model [11, IG, 30]. 

In general the toy model predicts significantly differ- 
ent /?aa vs. /aa curves as a function of centrality as a 
function of trigger momentum: the larger p^ (or, equiva- 
lently, smaller Kaway), the more concave the /?aa vs. /aa 
curve becomes. The concavity of the correlation is also 
in general larger for a CGC medium than for a Glauber 
medium. It would be interesting to see these predictions 
compared with future measurements performed over the 
full centrality range in small centrality bins. 

The strong suppression of the away-side jet should also 
lead to an anisotropy of the /aa relative to the RP. Fig. 5 
(a) shows the predicted correlation between /aa and 
^K. The corresponding correla- 
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assuming k, 



away 



,,-fAA 



tion between /aa and reduced anisotropy Wg'^* 
is shown in Fig. 5 (b). The results for Kaway = |k are 
shown in Fig. 5 (c)-(d). 

Several interesting features can be identified from the 
figure. First, v^^^ increases dramatically from P to P 
dependence, but even the P model under-predicts the 
PHENIX data [22] by about one standard deviation, 
where the uncertainty is currently very large. Second, 
calculated values are very sensitive to Kaway^ they 
are larger than the inclusive hadron V2 of Fig. 3 for 

Ka 



^away - but are less for Kaway 



jH. It will there- 



fore be very useful to see experimental results for 
as a function of trigger momentum. Finally, we also see 
a strong anti-correlation between i>2**/c and /aa, quite 
similar to that between U2/e and Rak- Such similar- 
ity may not be surprising if the physics of jet quench- 
ing for the trigger and away jets were identical (as in 
the present model calculation). A precision measurement 
of these two correlations, therefore, could either confirm 
such similarity or suggest new physics in the away-side 
jet quenching. 

The correlation between /aa and Uj"^* is also quite 
sensitive to the choice of the collision geometry: switch- 
ing from Glauber geometry to CGC geometry leads to 
about a 20% reduction of Wj** and vl^^ je at fixed /aa 
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FIG. 5: (Color online) (a) Centrality dependence of Jaa vs. 

from jet absorption model in 5% centrality steps and 
data (solid circles), (b) The same data and calculations, ex- 
cepts their ^2** have been divided by the eccentricities, (c), 
(d) are same as (a), (b) but are calculated for 



FIG. 4: (Color online) Centrality dependence of -Raa vs. 
7aa from the jet absorption model, compared with PHENIX 
data [17] in the 0-20% centrality bin. The calculations were 
done assuming k for away-side jet is (a) the same, (b) 1/2, 
(c) 1/3, or (d) 1/4 of that for inclusive jets. We compare the 
model curves in (b) to the left-most ([4-5] GeV/c trigger), (c) 
to the central ([7-9] GeV/c trigger), and (d) to the right-most 
([9-12] GeV/c trigger) PHENIX data point. The dotted di- 
agonal line indicates Jaa = Raa- The explicit predictions for 
0-20% centrality are shown as diamond and cross symbols. 

for Kaway = ^'^ (significantly smaller for Kaway = jk). 
Thus a precise measurement of this correlation may help 
in distinguishing between different initial geometries. 

Fig. 6 shows the correlation between V2 and ^2** for 
«^away = 5^ and Kaway = j^i- Sincc both obscrvablcs first 
increase then decrease from peripheral to central colli- 
sions, the correlation plot shows a rather sharp turn at 
around 20-30% centrality (iVpart ~ 150). This provides 
a rather precise way of identifying the centrality range 
at which the anisotropy reaches a maximum. The right 
panel shows the correlation of reduced anisotropics. We 
see that all four scenarios fall approximately on a univer- 
sal curve, especially for Kaway = 5^, but with different 
reaches along the curve. Apparently, the reach is larger 
for Glauber geometry and higher order / dependence. 
This implies that the efficiency with which jet quench- 
ing converts the eccentricity into anisotropy depends on 
collision geometry and I dependence. This curve also has 
an interesting shape: At small W2,r 0.2 for Kaway = 5^ 
and < 0.3 for Kaway = -jk), which corresponds to more 
peripheral collisions, V2 is larger than but then 

-Wj** > V2 at large V2,r (> 0.2 for Kaway = 5K or > 0.3 



for Kaway — 4K) ^. 

Before closing this section, we want to discuss all corre- 
lations together. In general, jet quenching models predict 
an anti-correlation between suppression and anisotropy 
(^AA vs. V2, and /a a vs. v^^^), while they predict a 
correlation between the suppressions (Raa vs. Jaa) £^nd 
between the anisotropics (^2 vs. W2**). Thus it is rather 
surprising to see from the data that Iaa has a larger 
anisotropy than that for i?AA (i'2** > yet it is less 
suppressed. This feature is very difficult to reproduce in 
our simple model (see Fig. 4 (d) and Fig. 5 (c)), and, in 
general, is a challenge for jet quenching theory [21]. 

Summary: The simultaneous description of multiple 
observables tightly constrains jet quenching models. In 
this paper we identified four correlations among single 
hadron and di-hadron observables as useful for constrain- 
ing jet quenching models: i?AA vs. V2, Raa vs. Iaa, Iaa 
vs. V2^^ and V2 vs. V2^^. Using a jet absorption model, 
we showed that these correlations are sensitive to various 
ingredients in the jet quenching calculations, such as the 
path length dependence, collision geometry, and input 
spectra shape. Specifically, i?AA vs. V2 is most sensitive 
to I dependence, i?AA vs. Jaa is sensitive to both spectra 
shape and path length dependence, and Iaa vs. V2^^ is 
sensitive to all three ingredients. 

We found that only our energy loss model with AE ~ 
P AdS/CFT-Hke energy loss and a CGC medium de- 



^ However for and CGC geometry, 112 > ?^2^* ™ centrality 
bins 
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FIG. 6: (Color online) (a) V2 vs. ^2*^^ ^-^d (b) i'2,r vs. «2** 
for 5% centrality steps and four cases for ft;away = ^f^- (c), 
(d) are same as (a), (b) but are calculated for 
The scaling seen in (b) and (d) is broken in more central 
collisions as the fluctuations in eccentricity e dominate the 

mean, (e> < J [e - {e)f . 



scribes the i?AA vs. V2 data well as a function of central- 
ity and is also qualitatively (within two to three standard 
deviations) consistent with the i?AA vs. /aa data as a 
function of the trigger px- While the P , pQCD-like en- 
ergy loss model is completely inconsistent with the i?AA 
vs. V2 correlation, the P model describes the i?AA vs. /aa 
correlations better than the P model, possibly indicating 
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